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1. INTRODUCTION 



The work reported here is part of a continuing program aimed at deter- 
mining the unsteady flow in a transonic compressor stage. The stage is 
installed in the Turbopropulsion Laboratories of the Department of Aeronautics, 
Naval Postgraduate School. 

This report has been compiled to facilitate use of the data acquisition 
and analysis programs which have been developed primarily for the study of 
unsteady fluctuating pressures on the casing inner wall. The equipment and 
instrumentation are discussed briefly in section 2. The data acquisition 
system and programs are outlined in section 3 and described in detail in the 
appendices. The post-real time data analysis programs are outlined with sample 
results in section 4 and are described in detail in the appendices. Section 
5 contains conclusions and recommendations for further work. 

2. EQUIPMENT AND INSTRUMENTATION 
2.1 The Transonic Compressor 

The transonic compressor test rig comprises an air turbine drive unit 
and an induction section which contains a filter, throttle, settling chamber 
and flow measuring nozzle. The turbine drive unit supplies 450 HP at 
30,000 RPM. The compressor is designed to operate at 30,460 RPM with a rela- 
tive tip Mach number of 1.5. At the design RPM and the tip Mach number, the 
flow angle is 65° and the pressure ratio is 1.6 at a referred flow rate of 
19 lbm/sec. The laboratory facilities and the test rig are described in 
detail by VAVRA and SHREEVE (1972) and VAVRA (1973). 
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2.2 Pressure Measurement 



Eight Kulite CQL-080-25 pressure transducers with natural frequency about 
125 kHz are mounted with their diaphragms flush with the inner case wall of 
the compressor. Further details are reported by PAIGE (1976). Table 1 in 
Appendix A gives the axial and circumferential location of the transducers 
relative to transducer number K6 which is the furthest upstream. The trans- 
ducers are used in conjunction with Datel Model 201C instrumentation 
amplifiers which have a flat frequency response to 100 kHz. 

Each Kulite pressure transducer is matched by a pneumatic static pressure 
tap at the same axial location in the case wall (except in one case - see Table 
1 in Appendix A) but displaced circumferentially. Other pneumatic static and 
total pressure taps are available upstream. A data recording system (VAVRA 
and SHREEVE, 1972) is used to record both the steady pressures from the 
pneumatic taps and the temperature data. The paper tape output from this 
system is processed using a Hewlett Packard Model HP9830A programmable calcu- 
lator to provide input data for the measurement of fluctuating pressures and 
to establish the compressor operating point. 

2.3 The Timing Disk 

To enable synchronization of the sampling of the pressure transducer 
outputs with the rotation of the rotor, an instrumented timing disk is fitted 
to the rotor shaft. The disk contains holes at intervals of one per rotor 
blade and one per rotor revolution. Light sensitive diodes and wave shaper 
circuits provide pulse trains to control sampling of the pressure transducers. 
This system is described in detail by WEST (1976) . 
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3. DATA ACQUISITION 



3.1 The System Hardware 

Figure 1 is a schematic of the data acquisition hardware with arrows 
indicating the flow of data and control signals. The system is under the 
control of the HP 21MX computer which operates either directly or through the 
device called "Pacer" to control the analog-to-digital (A/D) converter 
(model HP5610A) and which transfers data to the HP9867B mass memory unit via 
the HP9830A calculator. 




Figure 1. Schematic of the data acquisition system hardware 
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The peripheral device called the "pacer" is described in detail by 
WEST (1976) who originally called it RPACE. The pacer can trigger data 
acquisition from a stationary transducer at any fixed point in the rotating 
rotor frame, independent of the rotor speed. In effect, it divides the 
circumference of the rotor into 9 intervals, each with a circumferential 
length equal to that of the arc (measured at the wall) across two adjacent 
blade passages. Each of the 9 intervals are subdivided into 256 equal 
sub-intervals . 

The pacer receives one per revolution and one per blade input signals 
from the timing disk and performs two functions; it controls the timing for 
data acquisition and determines the speed of the rotor. 

3.2 The Program KULITE 

KULITE is the data acquisition program (in BASIC language) for the 
HP 2LMX computer. A flow diagram, listing, variable assignment and notes 
are presented in Appendix B. KULITE can be operated in three modes, viz. 
free- run, calibration and pacer. 

In free-run mode the A/D converter operates in mode 4 (see Reference 
No. 5). Up to 1616 samples of one A/D converter channel are taken with a 
frequency of 10^ samples per second. The sampling process is not synchronized 
with the rotor rotation. 

In calibration mode the A/D converter operates in mode 4 (see Reference 
No. 5). It scans through A/D converter channel numbers 1 to 12, taking 
1616 samples on each channel. The average of the 1616 samples is computed 
before the next channel is sampled. The scan is performed four times, with 
four different calibration pressures applied to the reference side of each 
Kulite transducer. 
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In pacer mode the A/D converter operates in mode 0 (see Reference No. 



5) . Sampling of a Kulite transducer output is synchronized with rotor 
rotation by means of two pulse trains generated from light beams chopped by 
the timing disk on the compressor shaft. One pulse train has a frequency of 
one per rotor revolution and the other has a frequency of 18 per rotor 
revolution; each pulse in the latter train corresponding to the passing of a 
blade past a fixed point. A full description is given by WEST (1976). 

In this mode a pressure transducer is sampled on successive revolutions 
at a fixed point in the rotating rotor frame. Currently, the sample interval 
is several revolutions of the rotor. Changes in program RPACE would allow 
samples to be taken at intervals of one revolution. If the flow can be 
regarded as steady in the rotating rotor frame this technique enables 
measurement of the wall pressure distribution "carried around" by the rotor. 
Flow unsteadiness in the rotating frame can be averaged or the frequency con- 
tent of the unsteadiness in successive samples can be examined. In this 
report only averaged data from 10 samples taken at each of 128 points across 
two rotor blade passages, is presented. 

In all three modes of operation the program KULITE transfers data from 
the HP21MX to the HP9830A. 

3.3 The Program TRAN4 

TRAN4 is the data acquisition program (in BASIC language) for the HP9830A 
programmable calculator. It receives data from the HP21MX computer, processes 
it and stores data on a disk of the HP9867B mass memory. A flow diagram, 
listing, variable assignment and notes are presented in Appendix C. 
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3.4 The Program RESET1 



RESET1 initializes a record number on the storage disk so that at the 
start of a run data can be stored in file DAT AY 1 beginning at the first record. 
The program is listed in Appendix D. 

4. DATA ANALYSIS 
4.1 The Program 

Off line data analysis is at present performed on the HP9830A with the 
BASIC language programs MAPI, MAP2, CONT, CONTI, PLOTSA, PLOTSB and TITIPK. 
These programs are described in detail, with listings, flow diagrams and 
notes, in Appendices E through J. 

MAPI is used to determine the sensitivity of the Kulite pressure trans- 
ducers from data acquired with KULITE in the calibration mode. In addition, 
MAPI is used to convert the voltages sampled at the pressure transducer out- 
puts to pressure coefficients. 

MAP2 is used to convert the 8 x 128 array of measured pressure coeffi- 
cients to a 29 x 128 array through quadratic interpolation in the axial 
direction. The program was written to reduce the effects of the course trans- 
ducer spacing in the axial direction. However, care must be exercised when 
it is used to interpolated across discontinuities such as shock waves and 
rotor blades. Linear interpolation is available through use of the program 
CONT or CONTI to plot contours of casing wall pressures. 

CONT is used to plot contours of constant casing wall pressure (in the 
frame of the rotor) from an array of pressure coefficients. (i.e. it pro- 
duces a wall pressure ,, map n ) . The program will accept any general rectangular 
array provided that the spacing in each direction is uniform. This latter 
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requirement restricts it f s use in this application to arrays obtained from MAP2. 

CONTI is used to plot contours of constant casing wall pressure when the 
array of measured pressure coefficients contains nonuniform spacing in the 
axial direction. Nonuniform spacing results in this application from the 
axial location of the Kulite pressure transducers. 

PLOTSA is used to plot (on the HP9862A plotter) the uncalibrated pressure 
distribution (in volts) across a blade pair for a given Kulite transducer. 

The input data is that originating from pacer mode of operation. The program 
is also used to plot the output of the one per blade signal from the timing 
disk. 

PLQTSB is used to plot (on the HP9862A plotter) the uncalibrated free- 
run data (in volts) from a given transducer against circumferential distance. 

TITIPK is used to superimpose the blade tip profiles on the wall 
pressure maps. 

4.2 Sample Results 

The results presented here are intended only to illustrate the capabili- 
ties of the programs. Comprehensive results will be given in a subsequent 
report. 

Figure 2 is a plot versus circumferential distance of the average 
pressure in the frame of the rotor across an arc of the casing wall equivalent 
to two blade passages. The pressure coefficient is defined in Appendix E. 

The plot was made with PLOTSA using data acquired in the pacer mode of operatior 
The precise location of each distribution relative to the rotor blades is not 
defined here. The locations are known approximately from the blade pair 
number specified in the acquisition program. They are located precisely from 
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Figure 2. Waveshapes of unsteady pressure distributions across two blade 
passages (uncalibrated and with arbitrary offsets) for the Kulite 
transducers. Data taken in pacer mode. 50% design speed; throttled 
to near surge. 8.7 lbm/sec referred flow rate. Pressure ratio = 
1.155:1. Blade pair #2. 
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-0.5 



- 1.0 




Figure 3. Typical record of raw one per blade signal taken in pacer 
mode. Blade passage frequency is 4.55 kHz. 




Figure 4. Typical record of data taken in free-run mode from Kulite number K10. 
The record is comprised of 1616 samples taken at a frequency of 
100 kHz. Each cycle is due to a blade passage with a blade passing 
frequency of 4.55 kHz. Compressor operating conditions as in figure 2. 
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Circumferential (opposite to direction of rotation) 




Figure 5. Contours of constant pressure coefficient ACp plotted 

by CONTI. (See Appendix E for definition of AC ). Blade pair number 
2. Compressor operating conditions as in figurl 2. 
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CIRCUMFERENTIAL 




STREAMWISE AXIAL 

Figure 6. Smoothed contours of constant pressure coefficient obtained from Figi 
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Circumferential (opposite to direction of rotation) 



-2-3 




Figure 7. Contours of constant pressure coefficient plotted by CONT with 
25x255 array. Blade pair number 2. Compressor operating 
conditions as in Figure 2. 
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Circumferential (opposite to direction of rotation) 
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Circumferential (opposite to direction of rotation) 




Figure 9. Contours of constant pressure coefficient plotted by CONT with 

8x128 array. Note there is a small error in blade location. Blade 
pair number 2. Sixty percent design speed. Throttled to near surge. 
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a knowledge of the orientation of the timing disk relative to the rotor in 
conjunction with the phase relationship between the one per blade timing 
signal and the pressure distributions from the transducers, A typical record 
of the one per blade signal, taken directly from the photo-diode output in 
the pacer mode is plotted in figure 3. 

The distributions shown are of output voltage from the Kulite transducers 
which at this stage have not been scaled to take account of the tranducer 
calibrations. The rapid changes in the signals, in particular those from 
transducers K9, KlO and Kll, are due to the passage of a rotor blade across a 
transducer and provide one means of estimating blade location. The distribu- 
tions have been plotted by linear interpolation between the 128 points across 
the blade pair. Each of the 128 points has been obtained from the average of 
ten samples, one sample being taken approximately each tenth revolution of 
the rotor. 

Figure 4 is a typical record of data taken in the free-run mode from a 
particular transducer. The fundamental frequency is that of the blade 
passage past the transducer (4.55kHz). The plot is a linear interpolation 
between 1616 samples taken at a frequency of 100 kHz. 

Figure 5 is a representative map of the contours of constant pressure 
coefficient with respect to an upstream reference static pneumatic pressure. 

It was generated with the program CONTI from the 8 x 128 array of measurements 
which were taken in the pacer mode and calibrated with the program MAPI. The 
map across two adjacent rotor blade passages is thus that of the mean wall 
pressures "carried around" by the rotor. The method of calculating the 
the pressure coefficients from the Kulite transducer outputs and the measure- 
ments from the pneumatic taps is given in Appendix E. 
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The location of the blade tip in figure 5 was determined from the circum- 
ferential location of the rapid change in pressure distribution associated 
with transducer K9 (figure 2). The transducer K9 has an axial location 
such that it is crossed by the leading edge of rotor blades. A more precise 
location of blades can be obtained from measurement of the circumferential 
location of the one per blade raw signal (figure 3) relative to the pressure 
distributions (figure 2) and from knowledge of the location of the one per 
blade holes in the timing disk relative to the rotor. 

Smoothing of the Wall Pressure Maps 

It is clear in figure 5 that linear interpolation between pressure 
coefficients in the axial direction across a blade gives incorrect contours 
along the blade between the transducer locations. This difficulty can only 
be overcome satisfactorily by using more Kulite transducers to provide a 
finer mesh than the rather coarse one provided by eight transducers. In this 
study the contours in the vicinity of the blades have been smoothed graphically 
by connecting with smooth curves those points in the map at which the pressure 
coefficient is the same and which lie on the lines scanned by the transducers. 
At operating conditions which give rise to shock waves it is possible that a 
similar smoothing procedure will be required. In the long run there is a need 
for a numerical interpolation technique which avoids interpolation across 
blades . 

Figure 7 is a wall pressure map produced by CONT for similar operating 
conditions to those used to obtain figure 5. However, pressures were 
measured at 255 points in the circumferential direction and MAP2 was used to 
generate pressure coefficients at 25 equally spaced stations in the axial 
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direction between transducers K6 and K12. It is clear that quadratic inter- 
polation does not avoid the problems of contour distortion by interpolation 
across a blade. In view of the very slow execution of CONT with a 25 x 255 
array of pressure coefficients, it is recommended that an 8 x 128 array be 
used with CONTI. 

Figure 8 is a wall pressure map of the same data as in figure 7 except 
that only every fourth circumferential point (of the original 255) was used 
in the reduction with CONTI (i.e. axial interpolation is linear through use 
of CONTI but MAP2 is not used). The resulting 7 x 64 array of pressure 
coefficients gives rise to very similar contours to those shown in figure 7 
for the finer mesh. In fact, the appearance of contours with = 0.4 near 
the blades in figure 7 might be an erroneous result of quadratic interpolation 
across blades. Note that the grid lines in figure 8 which indicate transducer 
locations are incorrectly plotted to have equal spacing, but this does not 
affect the above comments. 

Figure 9 is a plot of contours for a higher compressor rotational speed. 
The data reduction techniques were identical to those used to obtain figure 5. 
Note that the blade location shown was that calculated for the data of figure 
5. The location is slightly in error due to a difference in the tuning of the 
phase-lock loop circuit in the pacer for the second test. This problem has 
been solved recently by using the near-discontinuity in the pressure distribu- 
tion measured by the transducer at the leading edge to position the blade tip. 

Again it is stressed that the results presented are prelinimary and are 
intended merely to demonstrate the methods used and the capabilities of the 
system. 
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5. CONCLUSIONS AND RECOMMENDATIONS 



Programs for the acquisition and reduction of fluctuating casing wall 
pressures in a transonic compressor stage have been developed and run 
successfully with the compressor at this time operating at up to 60 percent 
of the design speed. Evaluation of the data acquisition system on a 
mechanical simulator indicates that it can operate over the full speed 
range of the compressor. In fact, at higher speeds the signal to noise 
ratio in the Kulite transducer outputs will improve significantly because of 
the higher pressures that will be encountered. 

The pacer mode of synchronized sampling has made it possible to deter- 
mine in a versatile manner the wall pressure maps in the rotating frame of 
the rotor. The pacer system can also be used to obtain measurements of 
flow properties away from the wall, e.g. flow velocity measurements with a 
dynamic probe. Wall pressure maps have been presented solely to demonstrate 
the capabilities of the pacer technique and of the data acquisition and 
analysis. Comprehensive data will be presented and interpreted in a 
subsequent report. 

There are some aspects of the programs which can be refined or which 
need further evaluation and the following recommendations are made. 

1. The subroutine RPACE in KULITE causes a sample to be taken in pacer mode 
about every tenth revolution of the rotor at 60 percent of design speed. 
This causes a delay in data acquisition which could be reduced by 
modifying subroutine RPACE. 

2. The degree of steadiness of pressure distributions in the frame of the 
rotor needs further investigation. This should begin with an examination 
of the standard deviations already computed in pacer mode. In separate 
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tests, a larger number of samples (at least 500) should be taken for each 
of several steps between blade pairs and the variations at each step 
examined for frequency content. The measurements should then be repeated 
with the case wall rotated peripherally by at least 90°. 

3. The technique of calibration of the Kulite pressure transducers under 
operating conditions effectively takes account of change in transducer 
sensitivity with temperature. Change in transducer offset (d.c. level) 
with temperature is presently handled by equating the time-average trans- 
ducer output voltage with the steady pressure obtained from a pneumatic 
static tap at the same axial location. The relationship between the 
steady pressure indicated by the pneumatic tap and the time time-averaged 
pressure at the tap needs further investigation. 

4. Because transducer K10 and pneumatic tap S10 are not coincident axially 
it is necessary to interpolate between readings at SlO and Sll. The 
interpolation in MAPI is at present linear but its adequacy has not been 
fully evaluated. 

5. The large pressure gradients in the axial direction across rotor blades 
are not resolved well because of the limited number of transducers. Two 
additional transducers, located midway between K9 and K10 and K10 and 
Kll, would greatly alleviate this difficulty. Linear axial interpolation 
across blades, as in CONTI, is misleading and hand smoothing of contours 
near blades is presently necessary. Quadratic interpolation, as in MAP2, 
does not solve the problem. Extrapolation of data up to but not across 

a blade surface should be investigated. 
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6* Shock waves have not been encountered at the low operating speeds at 
which the present data was obtained. The accuracy of resolution of 
shock waves should be studied in the light of the above discussion re- 
garding large pressure gradients across blades. 

7. The blades have been located on the wall pressure maps in this report 
from knowledge of the point in the circumferential pressure distribution 
(indicated by transducer K9 at the blade leading edge) at which the 
circumferential pressure gradient is steepest. This technique is subject 
to an, as yet, undetermined uncertainty due to irregularities in the 
geometry from blade to blade. The alternative procedure, whereby 
blades are located by use of the phase relationship between the one per 
blade signal and the circumferential pressure distributions, also needs 
further evaluation. 

8. In its present form the data acquisition system requires frequent key- 
board entries by the operator. In principle the system can be fully 
automated by pre-entering all necessary data with DATA statements and by 
replacing INPUT statements by READ statements. Some WAIT statements in 
KULITE would be needed to allow TRAN4 to catch up to KULITE. 

9. The format of graphical outputs can be improved by using the plotter 
to add alphameric information. 

10. Two-way data transfer beteen the HP21MX and the HP9830A is feasible. 

This capability should be developed to enable use of the faster HP21MX 
for repetitious data reduction. 
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APPENDIX A 



Table 1 contains the axial and circumferential location of the Kulite 
pressure transducers and the axial location of the pneumatic static pressure 
taps . 



Kulite 

transducer 

number 


Pneumatic 
static tap 
numb er 


Axial distance 
downstream of K6 
(inches) 


Circumferential 
location relative 
to K6 in 

direction of rotation 










K6 


— 


0 


0° 


K7 


— 


0.50 


+ 10° 


K8 


— 


1.00 


0° 


K9 


— 


1.37 


+ 10° 


K10 


— 


1.75 


0° 


Kll 


— 


2.12 


+ 10° 


K12 


— 


2.50 


0° 


K13 


— 


3.00 


+ 10° 




S6 


0 







S7 


0.50 







S8 


1.00 







S9 


1.37 







S10 


1.55 







Sll 


2.12 







S12 


2.50 


— 




S13 


3.00 


— 



Table 1 Location of Kulite pressure transducers and pneumatic static pressure 



taps. 
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APPENDIX B 



DETAILS OF KULITE 



KULITE is the data acquisition program for the HP21MX computer. Its three 
modes of operation are indicated in section 3.2. Figure 10 is a flow diagram 
of the program and a listing is given in Table 2. 



Variable Assignment for KULITE 



11 - Run //. Same as Run # in Log Book assigned to each 

start-up of the compressor 

12 - Test #. Refers to a particular operating condition 

within a run. 



13 - Day 

1 4 - Month 

15 - Year 

16 - A/D converter mode //. 

1 7 - Samples /channel in free- run mode. 

18 - Not used. 

19 - Experiment #. Refers to either 

(i) One time series of free-run data, 

(ii) A complete set of calibration readings (averaged) 
for all transducers, or 

(iii) Averaged pressures across one blade pair in Pacer 
mode . 

A1 - Channel //. Refers to A/D converter. 

T1 - Transducer //. 

N1 - Samples /point in Pacer mode. 

N2 - Blade pair #. 

M - Mean of pressure samples at a point in Pacer mode. 

S - Standard deviation of pressure samples at a point in 

Pacer mode. 

R - Mean of one/blade signal at a point in Pacer mode. 

L - Row number in K matrix of calibrations. 

A3,A4,A5,A6 - Associated with subroutine RPACE and defined by WEST (1976) 

A[101, 16] - Consecutive free-run samples stored row by row. 

B[101, 16] - Buffer in subroutine R5610 

C[10], D[10] - Buffers in subroutine R5610 
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E[4, 255] 



- - step 2 - - 



M. 



1 



M. 



3 



M, 



255 



S 



1 



S 



3 



S 



255 



A4 A4 3 



A4 



255 




R, 



‘255 



K[5, 12] 



This is the matrix of averaged data taken in RPACE 
across two blade passages. Note that 128 points are 
taken across two blade passages. This number can be 
changed to 255 by changing Line 195. 

Matrix of calibrations. Rows 1, 2, 3, each contain 
averaged calibration voltages for the twelve trans- 
ducers. Each of rows 1, 2, 3 corresponds to a 
different calibration level. Row 5 contains II, 12, 

13, 14, 15, 19 and the three reference pressures (which 
are keyed in on request) in K[5,7], K [ 5 , 8 ] , K [5 , 9 ] • Row 
4 is treated as another calibration level and is used 
to scan the offsets if needed. In that case any value 
can be input to K[5,10] for P Ref. 
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PACER 



Figure 10. Flow diagram for KULITE 



Cen& ) 
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Table 2. Listing of KULITE 



I F.EM "KULITE" DATALOGGING PROGRAM--SIMMONS-- 1 1 APRIL, 1977 
5 DIM AC 101, 1 6) , BC 101, 1 61 , CC 103, DC 101, EC 4, 2551, KC 5, 12] 

9 PRINT "ENTER RUN#, MONTH, DAY, YEAR" 

10 FOR 1=1 TO 5 

II FOR J= 1 TO 12 

12 LET KC I, J1 = 0 

13 NEXT J 

14 NEXT I 

15 INPUT 11,14, 13, 15 

16 PRINT "ENTER 0 IF 9830 IS RESET 1 FOR END" 

17 INPUT R8 

18 IF R8> . 5 THEN 47 0 

20 PRINT "ENTER TEST#, EXPERIMENT#" 

25 INPUT 12,19 

26 PRINT "0 FOR CALIBRATION OTHERWISE 1" 

27 INPUT Z 1 

28 PRINT# 8,*Z1 

29 IF Z 1< . 5 THEN 129 

30 PP I NT "ENTER MODE#, CHANNEL#, TRANSDUCER#" 

35 INPUT I 6, A 1 , T 1 

40 PRINT# 8; I 1 
45 PRINT# 8; 12 
50 PRINT# 8; 13 
55 PRINT# 8J 14 
60 PRINT# 8,‘ 15 
65 PRINT# 8; I 6 
70 PRINT# 8,* 19 
75 PRINT# 8,*T1 
80 PRINT# 8,’A1 
85 IF 16=0 THEN 165 
90 IF 16=4 THEN 103 
95 PRINT "WRONG MODE#" 

98 GOTO 3 0 

103 PRINT "ENTER SAMPLES / CHANNEL (NOT > 1616)" 

10 5 INPUT 17 
110 PRINT# 8; 17 

115 R56 1 0 (7, AC 1 , 1 1, I 7, A 1, I 6, EC 1, 1 1 ) 

116 PP I NT "DATA TAKEN" 

120 FOP. J = 1 TO 10 1 

124 FOR 1=1 TO 16 

125 PRINT# 8,* AC J, II 

126 NEXT I 

127 NEXT J 

128 GOTO 16 

129 LET L= 0 

130 LET L = L + 1 

131 PRINT "ENTER P REFCIN* H20 REL ATMOS) IF READY FOR CALIBRATION 
13 2 INPUT KC 5, 6+L 1 
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133 

134 

135 

136 

137 

138 

139 

140 

14 1 

142 

143 

144 

145 

146 

147 

148 

149 

150 

15 1 

152 

153 

154 

155 

156 

157 

153 

159 

160 

161 

162 

163 

164 

165 

166 

167 

168 

170 

180 

18 1 

182 

19 0 



Table 2. Cont. 



FOR J 1= 1 TO 12 

PRINT "STARTING CALIBRATION OF A/D CHANNEL # "J 1 
LET A 1 = J 1 

R5 6 1 0 <7, AC 1, 1 3 , 1 6 1 6* A1 * 4, BC 1, 13) 

LET KCL* J 13=0 

FOR J= 1 TO 101 

FOR 1=1 TO 15 STEP 2 

IF ABS(AC J, II )< .98 THEN 145 

PRINT "****•*■»****•»****■**■»********** " 

PRINT "****** OVERLOAD ON "J 1 "****** - 
PRINT "*■**■**«•*****■*■»■**■*•*■»**•»•»•****** ” 

GOTO 131 

LET KC L* J 1 3 =KC L* J 1 3 + AC J* II 
NEXT I 
NEXT J 

LET KCLjJ 13 = KCL.»J1 3/1616 
NEXT J1 

PRINT "0 FOR ANOTHER CAL I BRATI ON- -OTHERV I SE 1" 
INPUT Z 2 

IF Z2< . 5 THEN 130 
LET KC 5, 13 = 11 
LET KC 5, 23 = 12 
LET KC 5, 33 = 13 
LET KC 5^ 43 = 14 
LET KC 5, 53=1 5 
LET KC 5, 63=19 
FOR 1=1 TO 5 
FOR J= 1 TO 12 
PRINT# 8 J KC I j J 3 
NEXT J 
NEXT I 
GOTO 16 

PRINT "ENTER BLADE PAIR #* S AMPL ES/PO INT" 

INPUT M2* N 1 
PRINT# 8IN1 
PRINT# 3 * N 2 
LET A3= 0 

IF Tl>INT(Tl/2)*2+. 1 THEM 190 
LET A 6= 32768+ N 2*256-64 
GOTO 195 

LET A6=327 68 + N2*256 
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195 

20 0 

21 5 

22 5 

23 0 

240 

241 

242 

243 

244 

245 

246 

247 

25 0 

25 5 

260 

27 0 

280 

290 

300 

310 

320 

33 0 

34 0 

35 0 

360 

37 0 

38 0 

38 5 

39 0 

40 0 

41 0 

420 

430 

440 

450 

47 0 



Table 2. Cont. 



FOR 1=1 TO 255 STEP 2 
LET A3 = A6+ I 
L ET R= 0 

RPACE( A3; A4, A5) 

R5 61 0 C7, CC 1 3 , N 1 , A 1 , 0* DC 13) 

FOR J= 1 TO W 1 
LET 3C J, l 3 = CCJ3 
NEXT J 

LET A3=32768+N 2*256+1 
RPACEC A3, A4, A5) 

R5 610(7, CC 13, Nl, 0,0, DC 13) 

FOR J=1 TO Nl 
LET R=R+CC J 3 
NEXT J 
LET R=R/N 1 
LET S= 0 
LET M= 0 
FOR J= 1 TO Nl 
LET M=M+BC J, 1 3 
NEXT J 
LET M=M/N 1 
FOR J= 1 TO Ml 

LET S=S+( (EC J, 1 3 -M>* ( 3C J, 1 3 -M) ) 
NEXT J 

LET S=SQR(S/(N1- 1 ) ) 

LET EC 1 , I 3=M 
LET EC 2, I 3 = S 
LET EC 3, I 3 = A4 
LET EC 4, I 3 = F. 

NEXT I 

FOR J= 1 TO 4 

FOR 1=1 TO 255 STEP 2 

PRINT# 8,*ECJ,I3 

NEXT I 

NEXT J 

GOTO 16 

END 
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Notes on KULITE 



1. A/D converter channels . It is essential that the "raw" one per blade 
signal be input to channel 1 of the A/D converter. Allocation of the 
other channels is not unique but the allocation in Table 3 is recommended. 
Channels 11 and 12 are scanned but at present are not used in subsequent 
analysis. 



A/D Converter 
Channel Number 


Signal 


0 


one per blade raw signal 


1 


K6 


2 


K7 


3 


K8 


4 


K9 


5 


K10 


6 


Kll 


7 


K12 


8 


K13 

* 


9 


P , - P 
ret atmod 


10 


S2 - P . 
ref 


11 


Unused 


12 


Unused 



Table 3. Allocation of signals to A/D converter channels. 

* 

P -is pressure applied to reference side of Kulite transducers, 
ret 



2. Subroutine R5610 is described by WEST (1976, p. 17). 

3. In calibration mode the scan through the twelve channels must be made four 
times. The first scan must be with- the pressure tapping S2 applied 
simultaneously to the reference side of Kulite transducers. The second 
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and third scans must be made with other steady pressures applied to the 



reference side of the Kulites 0 The fourth scan must be made to satisfy 
the program but at this stage the data taken is not used in subsequent 
analysis. This scan is included to enable logging of the offsets on the 
Kulite amplifiers should they be of interest. 

4. The program searches for overloads (i.e. greater than 0.98 volts or less 
than - 0.98 volts in the calibration signals). If it detects an overload 
among alternate samples in the 1616 samples taken from any transducer the 
offending A/D converter channel number is displayed, the scan is aborted 
and the program is reset to repeat the scan. The limit of 0.98 volts can 
be changed in line 140.. 

5. The one per revolution signal from the timing wheel indicates the origin 
for circumferential measurements around the rotor. The pacer then uses 
the one per blade signal to divide the rotor circumference into 9 equal 
intervals, the first interval beginning at the origin. These intervals 
are designated by blade pair numbers, although the start of an interval 
need not coincide with a blade tip. Each interval represents a circum- 
ferential length, in the rotor frame, equal to that of the arc (measured 
at the wall) across two adjacent blade passages. Each of the nine 
intervals is divided into 256 sub-intervals. In pacer mode the scan 
across an interval begins after the first sub- interval and ends after the 
255th sub-interval. With stepping sequentially across the sub-intervals 
in pairs, a total of 128 points are sampled. It is convenient to take 10 
samples at each point (one sample approximately each ten revolutions) to 
compute the mean and standard deviation. 
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6. Even numbered transducers (e.g. K6, K8 etc.) are located on one axial 
line and odd numbered transducers are on another axial line which is dis- 
placed around the casing wall by 10 degrees in the direction opposite 
that of rotation of the rotor. The parity of the transducer number is 
evaluated in line 180. 

The variable A3 determines the time (in terms of degrees of rotation 
of the rotor) after the one per revolution pulse when a sample is taken 
at point I. For example, in line 190, A6 = 32768 4- N2 * 256 for odd 
transducer numbers, and A3 = A6 4- I . This defines the sampling time 
(approximately each tenth revolution) for point I in blade pair N2. 

(I = 1 to 255, N2 = 1 to 9) . Point I can be sampled 10 degrees earlier 
for even numbered transducers by setting 
A6 = 32768 + N2 * 256 - 64 
and A3 = A6 4- I 

The subroutine RPACE is described in more detail by WEST (1976) . 

7. During each scan of a transducer across a blade pair the raw one per blade 
is also sampled. This is used later in TRAN4 to determine the location 

of the measured pressure distribution relative to the rotor. 
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APPENDIX C 



DETAILS OF TRAN4 



TRAN4 is the data acquisition program for the HP 9830A programmable calculator. 
Figure 11 is a flow diagram and a listing is in Table 4. 



Variable Assignment for TRAN4 



N1 

Z1 

A3 

A2 

T1 

T2 

A[5,128] 



RECORD it. i.e. Number of first available record in DATAY1 

- IDENTIFIER ( = 0 FOR CALIBRATION - OTHERWISE 1) 

- BLADE PAIR it 

- SAMPLES /POINT in Pacer mode. 

- LOCATION (between 1 and 255) of point in blade pair where 
1/blade signal is 0.5 volts and increasing. 



Location of point in blade pair where 1/blade signal is 0.5 
volts and decreasing. 




M. 



step 2 



M 



255 



S 1 S 3 S 255 

\ 1 A 4 ^ A4 255 

R 1 R 2 R 255 

II, 12, 13, 14, 15, 16, 17, Tl, Al, A2, A3, 0, 0, 0 



NOTE that variable names in this array are those used in 21MX 
program. They should not be confused with variables used in 
TRAN4. This is the matrix of averaged data taken in RPACE 
across two blade passages. Note that 128 points can be changed 
by changing dimension statement and the FOR loop. 

B [ 9 ] - Buffer for identification data 

[II, 12, 13, 14, 15, 16, 19, Tl, A1] T 

NOTE that the variable names in this array are those used in 
21MX program. 
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D[l6> 102] 



K[5, 12] 

RECY # 
DATAY1 



Free run data 
Stored row by row 



I II 

I 12 

1 13 
*14 
I 15 

lie 

1 19 
* T1 

| A1 

I 17 

Column 



I 



102 



NOTE that variable names in column 102 are those used in 
21MX program. 

Same as K[5, 12] in 21MX program. 



DATA FILES 



This file contains 1 record. It contains N1 which is the 
number of the 1st available record in file DATAY1 

Data file for calibration, free run and Pacer data. Le. for 
K, D and A arrays respectively. It contains 300 records. 
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PACSR 
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Table 4. Listing of TRAN4 



l :L.I ' I Klin i • .. [ IT!. 'No * -• !. 7 i : i !■ ' 



1 ! ! | • 


i 

. L.h'.' 


h >; 


- i.i m r at i 


48 [j 


i i‘i | i 


si: i 4 i j 


m 2 1 ? h Of. 5 ■ 1 


!: & R 


ER.0 


l! 1 ' N 1 




bU M 


i :: i r i) 


- 2ER 




78 MOT A 


~2ER 




80 R 


Ell T 


EOT IF 


NEXT jOT A 


9 0 E 


N TER 


* i 'i -f .> 


•l 

'. i. 


1 0 0 


1 1- 2 


1 . 0 „ 8 


THEN 2O0 


i '!. 8 


PEM 


RFCE 5 V 


E I BENT I F I 


1 20 


FOR 


1 = 1 TO 


9 


1 30 


ENTE 


f. 1 j ■* 


> 6C I 1 


140 


NEXT 


I 




180 


R E M 


TEST M 


ODE# 


1 6 0 


IF B 


[ ] = 0 


THEN 530 


1 “7 m 
i i r. .' 


IF B 


C 6 J = 4 


THEN 360 



180 HI SR "WRONG MODE #"5 

1 " 8 R E M R E C E I V E 8 N D S I 0 R E C A l . i B R A 71 0 r l 

200 FOR 1=1 TO 5 

210 FOR J=1 TO 12 

2 _ O L N T E K 1 •> * ■ K l. I j J ] 

230 NEXT J 

240 NEXT I 
250 PRINT 

260 PRINT "HRRRY OF CALIBRATION VOL TAG 

270 PRINT 

280 NOT PRINT K 

200 RE OB # 2 s N 1 

380 NOT PRINT # 25 R 

310 Mt=Hl+l 

326 READ HI 5 l 

38 PRINT HI? Ml 

340 GOTO 920 

350 REM RECEIVE BHD STORE FREE-RUN LOT 
360 ENTER Cl f *>17 
370 FOR 1=1 TO 16 
380 FOR J=1 TO 101 
3 '90 ENTER 1 1 > * BC I j J ] 

480 NEXT J 

410 NEXT I 

420 FOR 1=1 TO 9 

430 DC I j 102 ]=BC I 3 

448 NEXT I 

458 DC 10! 102 .1=17 

4 6 U R EBB # 2 s N 1 

470 MAT PRINT * 25D 

438 Nl=Hlrl3 

490 READ #1»1 

50Q PRINT #15 N1 

510 GOTO 350 

520 REM RECEIVE PACER DATA 

530 ENTER Cl, OA2 

5 4 0 ENTER C 1 * ■' A 

550 FOR J=1 TO 4 

560 FOR 1=1 TO 123 

5-f, ENTER Cl,r,flC J, I 3 

580 NEXT I 

590 NEXT .J 

60O REM FIND POSITION OF REFERENCE FIJI. 



[i 

; 1 



BfiT n 



ES ’ 



fi 



SE •]. ELhDE 
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Table 4. Cont 



h. : :■! 


I - : U 




6 £:• 


[ - 1 + 1 




l;i 9 0 


IF i ' 1 


Hfc'il i i- 


6 4 0 


j. 1 H L 4 ■' !. ! 


-o." HLn ... 


550 


GOTO o20 




6 G 0 


T 1 ■■■■ 1 




670 


1=1 + 1 




680 


IF AL4.11 


• . 8 , : ' f i ll T 1 ; : 'i 


690 


g u r u g 7 y 




7 0 0 


T 2 1 




710 


•. 3 U 1 U i ■ Pj 




?;::0 


I) ISP • REFERENCE PU! Si 


730 


FOR 1=1 r 


0 9 


740 


AC 5 » I J— BC 


I I 


750 


NEXT 1 




760 


1 1 [ j !» 1 y ] “ r i 




770 


AC 5 . 1 i J -A 


3 


7 8 0 


9 EM 8 * •Ji r ’E 


P M G E F DH Th 


790 


READ #2 .. N 


1 


8 0 0 


MAT PRINT 


# 2 S H 


8 1 0 


N 1 = fi 1 + 5 




820 


READ # 1 ? 1 




8 3 0 


pp irir tti; 


Nl 


o 4 U 


REM P Flil T 


E X )~ L F I i i E N f 


850 


PRINT 




y i“i y 


PRINT 




879 


WRITE a? 


i 8 y B L I ] ? B 1 


830 


FORMAT S'! 


* “RUlUr *F5 j 


690 


HR INI ' MUDi.: ti : BL o ] !» 


900 


IF BC 5 1 = 8 


THEM 1 050 


9 1 0 


IF 2 r- 0 . 5 


THE M 9 9 0 


920 


PRINT 




9 3 u 


PR I NT 




9 4 0 


WRITE 15 


■i UE'U 1 E i.. 3 ’ 1 J 


9 5 0 


PRINT "EX 


PEP i WENT# '■ Rl 


960 


PRINT "Cfi 


L I BRAT I C N — 


970 


PRINT " 




9 8 0 


G 0 T 0 1 8 1 8 




990 


PRINT' “TR 


A NS DUi-ER 


188 


0 PRINT "S 


AMPLER, C HAH 


1010 PRINT 




1 82 


0 PRINT 




183 


8 PRINT 


-J- t -T- A -A 4- Ar h +■ A -+• 



.1 i. i 



no f r ou Hlv 



EXPERIMENT #"eC 7'3 



" , F 



5 i> t j 



RECft " M 1 - 1 



U H h N N E L # ‘ 6 1.. 9 1 
EEC tTNl-13 



1U90 GOTO 1120 

1050 P R I H T ' T E: A h 0 DU CEE # 1 Bf 8 ] > “ C H A N N EL # " B C 9 I 

1 0 6 0 W R I T E 1 .1 5 ' 1 0 8 0 > fi 3 > R2 j H 1 - 5 

1070 PRINT " REFERENCE PULSE AT POINT'#" FT 1+T2 ■' 2 
1 080 FORMAT " BLADE PA I R * “ > F4 . 0 » 5X j " SAMPL.ES/PO I NT " , F4 . 8 , 
1090 PRINT 
1100 PRINT 

1 i 1 t-T F' R' I [ \ I ' 1 -r- -4 •+• "l- “(r •+• -A ■+• -.i- 7 A 7 -* ■+• +• -A 77 ,T 7*7 *?:' 7*7 +• 7*7 7'7 7*7 77 -V ’+■ +■ p 7*,- ■+• T '#7 ~*r "+* 

1120 D I SP "ENTE R 1 FOR MORE DATA" 5 

1130 INPUT Q 

1140 IF 0=1 THEN SO 

1150 PRINT 

1150 PRINT 

1170 PRINT 

1180 END 
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Notes on TRAN 4 



1 . 



2 . 



(line 30) DATAY1 is a file on a removable disk for temporary storage. At 
the end of Run number n (n is a two digit integer) a file CKRWn must 
be opened on the fixed disk and data must be copied into it from DAYA1 
for long term storage. The number of records in CKRWn must equal the sum 
of record number printed out with identification of the last experiment 
and a number k where 

k = 1 if last experiment was a calibration 

=5 if last experiment was in pacer mode 

= 13 if last experiment was in free- run mode 

Figure 12 contains sample print-out from TRAN4 for the three modes. 

(Line 600) The position of the centre of the raw one per blade pulse is 
found relative to the pressure distribution across a blade pair by 
searching through the 128 averaged samples at each point to find the 
sample which first exceeds 0.5 volts and the first subsequent sample 
which is less than 0.5 volts. The corresponding point numbers are 
averaged. If the pulse is not found (due to inadequate signal level) , 
the program displays an ERROR. After the correct one per blade signal 
has been re-established both KULITE and TRAN4 must be rerun. 
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Figure 12. Sample print-out from TRAN4 for calibration, mode 0 (pacer) 
and mode 4 (free-run) operation. 



i\iji - ■ 



E Rill limit I 
i BF'OT 1 ON - 



! L • . I Jr 



PC C It ' 






4 j; -r- “I" -4. V +• k h V 717 *- -k ;h '-Jr -h 4 -i • 4 -i; 7; H I- •! »• -!;• •+ .... I- V >. •+• ~k 717 -t~ jk •;,* 4 4; 4| 



j-’Ur i ri 58 TEST# i 

r • It !.■ L # 4 EXPEL i MEN! # 

I i • ! t“i !' J Di.l’/EP rt o 



CHFiHNEl tt I 
REG tt 2 



k 4“ ■+■ -r- +• 7 i- -4" -i" 717 7k •+" ;i“ 1 -i I7 717 ;r •+• "17 -r 17 717 747 *;i; •+• -*7 4* 47 4- **!•• •+• -f; • t- -f • «r 



f ; : . : D t : ,t Q 

1 r’.ri!’- * iJiJL L_ ! • ' 4 »: ■ 
L* L n L.i It., r t' 1 i K t* 

REFER E ML E 'RJLR 



TEST# 1 5- 24' 1977. 

£ X r t R [ M £ H I # -8 

CHFiNHEL # 1 

8 AMPLE S--'' P U I I t 1 1 0 k L t # i 5 

HT POINT# 31 



•-r* -»• -n- •+• 4 -r ;i- +■ ;i; -j;- t- 717 4 i • • •+ ;• 717 -77 •+■ -4* -k -17 77 ■+■ -#7 #7 4“ ■+■ 4- 77 4- 717 717 77 717 -4- 7(7 4*> •*(— ^ 717 717 •+• #7 7(7 -4- -4»* 4- -+• 717 -<7 1- 4« • 47 •••;• 7I7 t" 



P u 1 1 # 


j C 


TEST# 1 




NO EE # 0 




EXPERIMENT 


tt 4 


T RAN ’ DU CEP 


If b 






B L. h D 1 . P H 1 F 


4 


i SAMPLES/ Pi 


01 NT 


REFERENCE ; 


rULXt 


HT POINT# 31 




+• -1- -77 r* 77 "17 -k +* 7!" • k 






+■ 4k -k 77 ■ 


PUN# 




TEST# 1 




MODE # 4 




EXPERIMENT 


# Si 


i R ! i N 0 D IJ L E k! 


# ? 






SHHPL.ES/Ch 


El Nth EL 


1 6 1 6 





CHFIHNEL # 1 
10 PEC tt 20 



717 7)7 -77 717 “k -I- 7k 717 4- 7I7 -+■ 7I7 4- 4- 47 ■+• 7I7 7I7 •+' 717 4” 77 •+■ -jk 717 -i- 4- •+• -V -4- -I- 4, 



5 24 1977 

CHANNEL # 2 
PEC # 25 
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APPENDIX D 



DETAILS OF RESET1 



RESET1 initializes to 1 the number in the file named RECY // . This 
enables data acquired at the start of a run to be stored at the start of the 
file named DATAY1 . Table 5 is a listing. 



Table 5. Listing of RESET1 



?. \< L ' ■ ! "RE S Eli" - S I M M 0 N S 8 E T 3 R E C V* T 0 l 

E i ::s uiilT' u 

30 F IlDj REC Yfe 

4 >-1 H 1 

5 0 R £ H .[I # 1 .i 1 

i- O PRINT # 1 ? d 
, ' U kb R i.i # 1 j 1 

£0 READ # 1 5 H 
90 PRINT H 
1 0 0 £ \ \ D 
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APPENDIX E 



DETAILS OF MAPI 

MAPI is used to compute the sensitivities of the Kulite pressure trans- 
ducers and to convert voltages sampled at the pressure transducer outputs 
to pressure coefficients. Figure 13 is a flow diagram and a listing is in 
Table 6. 

Variable Assignment for MAPI 



N 


compressor RPM 


T 


total temperature (called elsewhere) measured at axial 

location of S2 (entered in degrees F) . 


PI 


Static pressure measured at S2 (entered in inches of 

water absolute) . 


P2 


Total pressure (V measured at axial location of S2 

(entered in inches of water absolute) . 


M 


square of Mach number at axial location of S2. 


U1 


2 2 

square of rotor tip speed in ft /sec . 


A1 


square of speed of sound at axial location of S2 in ft/sec. 


Q 


Reference pressure for computing pressure coefficients. 


Cl, C2 , C3 ■ 


First, second and third calibration pressures applied to 
reference side of Kulite transducers (inches of water 
relative to atmospheric) . 


RO 


Record number for calibration data on disk. 


TO 


- Kulite transducer number. 


A 


Kulite transducer sensitivity in inches of water per volt. 


PO 


Pressure from pneumatic wall static tap corresponding in axial 
location to a Kulite transducer. (inches of water relative 
to pressure at S2) . 
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SET QoARTeR — 


<P 


— Tt 
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Figure 13. Flow diagram for MAP 1 
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Table 6 



Listing of MAPI 



.... k .': i k : .i i 1 

3 }j I-* !, l ;i '-'T ! 



7k* h 88 J ui I fcl-t- - 1 8 
0 hi jl i : . p ' jf I' | f !: : <•••' f • : i 
9 ! "' IftbUT r, 



iijk 


oi;>p 


' EM if.' 


■ i ; !■ ■ , n* ; 






1. 1 kj 


!. i ! P 1 1 T 


r 








120 


j - [ h-J 


: : N 








130 


: i i. '3 


"ENTER R 


V t i , f | ifi.il i „ \ i , ■; r. 


hh •: " 3 




i ; 0 


i i if ij ! 


i • i . r :■ 








100 


,'i -* 1 < - : 


; ■ p j > !•.. H 


••7 !. , -1 1 - L :« ▼. 3 -'0 4 






100 


. j i _ R 


i ‘N-.N:.. 5 


. 6 U ! * 2 






1 "0 


r • J - 1 . 


5 >. ;4- 


i -N . 3 - i ■ 1 >-02 . 






180 


fji - 0 _ / 


- F U .IMI i 


1 ill. ' 






1 80 


D 1 3 f 


" I: ri T L ! ■: 


L rl L PRESSURES - 


I f t .. h 2 0 K: E 


L 1 1 f i < 


8 i. ; j 0 


I riPU T 


Pi ? 03 « C 


3 






2 1 0 


:.i I : r .p 


"ENTER i 


4!. REl It 3 






220 


INPUT 


R 0 








2: 3 0 


READ 


4 J. > P 0 








24u 


MR T R 


END tt ! 5 








250 


r 0 •= 6 










266 


PRINT 


' TRhiTP. 


HO. - "T0 






£ 7 0 


r j i j 3 1 j ? 


' 4 !’ 








2 8 0 


PRINT 


* 3 E !• 3 . 


OF T " T 0 " = A "IN. 


H 2 Li VOLT 




2 8 


HI TP 


'EM TEE f 


Mi 'ILL FOR T"T8"[H 


ri 8 0 REL 


8 n 


*3 0 0 


INPUT 


P 0 








3 1 U 


PEN i 


P TO-- N 


[HTEPPuLflT r OH IS 


NEEDED E 


R 1.. . | i i J j::. 


3 2 vj 


IF HB 


3 i‘0- L U ■ 


>0.5 THEN 360 






3 36 


n <_ 3 p 


“Also eh 


ITER PH ALL FOR Tl 


i " 3 





340 INPUT P 9 



350 


F 0 ~ p U r 1. 


H 4 - 


f ' 0 ■ - U 


. 2- 9 


n I* 






360 


DISF "E 


NIL 


R PHC 


ER c 'l 


RE 


33 RE 


OItt FOP ! TO; 


3"-0 


INPUT R 


i 












3 80 


READ ffl 


i R 1 












3 9 0 


MAT PEA 


D tt 


I •: A 










4 0 0 


NO =0 














4 i 0 


FOR I - i. 


T 0 


32 










420 


i — \ 
i — . 
i i 


Af 1 


• NO 4 3 


2-f I .3 


- c 


C 1 .. It 


j - 5 3 > H - Fl + P 0 > . C 


430 


NEXT i 














4 40 


READ #2 


i T 


rj 6 ) - 


4 -NO 


I 1 .1 






4 5 6 


MAT PR I 


NT 


# 25 P 










460 


L.[ TO- 5, 


NO + 


I 1=F'C 


1 J 








470 


FOR 1=2 


TO 


O ‘T 1 










480 


IF PC I 3 


>LL 


T 0-5 ? 


HO-*- 1 


3 


THEN 


500 


4 9 0 


LC TO -5 - 


NO + 


1 T=P[ 


I 1 








500 


NEXT I 














5 1 8 


!JC TO -5- 


NOi- 


1 J=PC 


1 3 








520 


FOR 1=2 


TQ 












530 


IF PC I 'J 


HE 


TO- 5 ? 


NU-l 


3 


THEN 


550 


549 


UC TO-5 < 


HO + 


1 ] = PC 


I 3 








550 


NEXT • I 














560 


rlO=NO J - 1 














570 


IF NO 3 


C ' 

. t 


THFN 


OHO 








580 


SOTO 4: 


0 












5 8 0 


TO=TOt 1 














6 0 0 


IF T 0 : .1 


j .. !..* 


THEN 


y\ 









1 y 'lii i l 

* ■•* U; ■■ ’• ll 

• k . «' r f' 8 1 



3 I U hPE Uh F 31: 



45 



Table 6 



Cont 



. ! { ) j i . v I 

» « l p =•; | ; 

; OK I ■ T. i 

"UK • ! • 1 TO 4 

i F HI 1 ? J 3 U ! HL'i I i- F"V! 

u-ur i , I j 

NEXT J 

NEXT 

PRIM'! "CP 

I -- i l ; : 0 

FOR i-1 TO 0 

FOR j--l TQ 4 

IF LL I :• J i " i- 3 Hi; hi >‘5U 

L r - L C i .« J 3 

NEXT 

NEXT I 

PRINT "CF i‘i 1 '■ L 

STOP 

PEN .uBk'Jij ; .1 ML 10 CRLuULhlE iRMNSDUUEP SEHSi i i l T Y 

8 - S' * ■ C J *• C [ ; 7 O •• ' 3 C 2 * C [ 2 .« T 0 5 3 + C 3 C C 3 > f O - 5 .] .> 

8 ~H " ' ' j r i _ L ;; i- I ; .• X *. ij !_ | i 7 hj — hj 1 r IJ [ A f hi — 0 J "I” !_ 7 i 3 ij — .j ,] ’ 

H - fi . < C 1 ' 2-'-02 l 2-+-C3'T'2 > - ' C 1 +C2 4-C3 > 3 2 > 

h - - 1 R 

RETURN 

END 



I hi . H , 
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are 



R1 - Record number for pacer data Q 

NO " NO+lis pacer data quarter number. The 128 samples 

divided into 4 sets of 32. 



U 



Maximum pressure coefficient in set across two blade passages. 



L 



Minimum pressure coefficient in set across two blade passages. 



AS [5, 128] 



Same as A[5, 128] in TRAN4 . 



C [5, 12] 



Same as K[5, 12] in TRAN4. 



P [32] 
U[8, 4] 



Array of pressure coefficients for one transducer and in one 
quarter of pressure distribution across two blade passages. 

U[I, J] is local maximum pressure coefficient for transducer 
number I and quarter J. 



L[8, 4] 



L[I, J] is local minimum pressure coefficient for transducer 
number I and quarter J. 



DATAY1 



Same data file as in TRAN4. 



PRESS - contains pressure coefficients. The 32 records contain in 

order quarters 1, 2, 3, 4 for transducer 6, quarters 1, 2, 
3, 4 for transducer 7, etc. 



Notes on MAPI 

1. The correct CKRWn file name for the run under consideration must be 
entered# 

2. Storage of pressure coefficients . The 128 pressure coefficients asso- 
ciated with each transducer have been grouped into 4 sets (quarters) with 
each set being stored in a separate record. This has been done to 
facilitate use of interpolation programs such as MAP2. Interpolation 
expands the size of the data array so that pressure coefficients must 

be recalled from the mass memory in subsets in order to meet the storage 
limitations of the HP9830A. 
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3. Sensitivities of the Kulite transducers are computed as follows: 



Cl, C2, C3 are the three steady pressures applied to the reference side 
of the Kulite transducers. They correspond to mean output voltages 
C [1 , TO -5], C[2, TO -5], C[3, TO -5] from Kulite transducer number 
TO. The sensitivity of a transducer (i.e. the slope of its calibration 
curve at a particular mean operating temperature and pressure) is ob- 
tained from a least squares fit of a straight line through the three 
points . 

4 . Calculation of pressure coefficients . 

The Kulite transducer output voltages E are converted to pressure 

coefficients C as follows: 

P 

C = ( (E - E) * A 4* PO) /Q 
P 

where E = transducer mean output voltage obtained during calibration 
with S2 on reference side of diaphragm. 

A = sensitivity of transducer in inches of water per volt. 

PO = mean wall pressure (at same axial location as Kulite trans- 
ducer) measured with pneumatic tap. (inches of water relative 
to S2) . 

Q = reference dynamic pressure (inches of water absolute) computed 
as follows: 

The reference dynamic pressure is expressed in terms of the upstream 
density p and the upstream flow velocity measured in the rotating frame 
of the rotor. Hence 



1 2 2 
Q = j P (v 4* u ) 

where V = flow velocity at station S2 (ft/sec). 
U = rotor tip speed 



ttN 5.5 ~ 

30 X 12 ft,sec ’ 

N - rotor RPM 



(i) 
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It follows (noting that variable names are not necessarily the same as 
in the listing of MAPI) that 

Q = \ YP (M 2 + U 2 /a 2 ) (2) 

where P, a, M are static pressure, speed of sound and Mach number at 
station S2. 




where t t and are total temperature and pressure at S2. 

By introducing (1), (3) and (4) into (2), Q can be calculated in terms 
of P, P T and T t . 

5. Note that S10 and K10 are not at the same axial location. For purposes 
of computing pressure coefficients an effective mean wall pressure at K10 
is obtained by linear interpolation using values at S10 and Sll. 

6. The maximum and minimum pressure coefficients are computed to aid in 
the choice of contours when using the programs CONT or CONTI. 
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APPENDIX F 



DETAILS OF MAP 2 

MAP2 is used to convert the 8 x 128 array of measured pressure coeffi- 
cients to a 29 x 128 array through quadratic interpolation in the axial 
direction. A listing is in Table 7. 



Variable Assignment for MAP 2 



B [ 32 ] 



temporary storage of pressure coefficients 



C[29,32] - array of interpolated pressure coefficients across one 

quarter of a blade pair. 

P[8, 32] - Array of pressure coefficients at Kulite transducer locations 

across one quarter of blade pair. 



X[8] 



Axial location of Kulite transducers downstream of trans- 
ducer K6. (inches). 



Q 



quarter number. 



PRESS 



Same data file as in MAPI. 



INTER 



File for storage of interpolated pressure coefficients across 
one quarter of blade pair (15 records). 



Notes on MAP2 



1 . 



Lagrangian interpolation is used, i.e. if P^, P^, and P^ 



are known at 



x 1 , x 2 , x^ then 



:■ limp 



Table 7 . Listing of MAP2 



— • M 



Oj 



P-- 







uJ 





Li_ 













JJ ! : ! 
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•I - 1 -■ 



T.il-lc 7. Coni. 
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b 4 W ML X ! ! 

b L :.n REN STORE } NT ERPOLRTLU PRESSURE COER" I f I F II f S Oh HI 
6 s 0 RE fill m j« i 

c?q Mm t pr nr # 2;c 

680 TUSP "END" 
bOO END 



2. As a result of the equispaced interpolation, the transducer axial loca- 
tions will in general not coincide with axial stations in the interpolated 

array. 

3. Three interpolations are made in the axial direction between transducer 
measurements . 

4. The program must be run for each quarter of the array of measured 
pressure coefficients. After each running of the program the contours 
must be plotted with CONT prior to running MAP2 for another quarter. 

5. Interpolation between X[l] and X[2] is made with a quadratic through 
pressure coefficients at X[l] , X[2] and X[3]. Interpolation between 
X[7] and X[8] is made with a quadratic through X [ 6 ] , X[7] and X[8] . 
Interpolation between X[I] and X[I + 1] (for I > 1 and < 6) is made by 
averaging the quadratic through X[I - 1], X[I] and X[I + 1] and the 
quadratic through X[I], X[I + 1] and X[ 1+2]. 
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APPENDIX G 



DETAILS OF CONTI 

CONTI is used to plot contours of constant casing wall pressure from 
an array of measured pressure coefficients. The program handles the non- 
uniform axial spacing of the Kulite transducers but requires uniform 
circumferential spacing in the array. CONTI is written to accept the array 
generated by MAPI. Figure 14 is a flow diagram and a listing is in Table 8. 

Variable Assignment for CONTI 

P[I, J] - array of pressure coefficients 

A[3, 3], B [3 , 3], Z [3] , Q[3], D[2, 2], F[2], G[2, 2], H [2] - defined in 





note 6 of this appendix. 


x[9] 


axial location (inches) of Kulite transducers downstream 
of K6 . 


10, JO 


dimensions of P[I, J] in axial and circumferential directic 
respectively. 


c 


value of pressure coefficient on a contour. 


E 


triangular element number. 


El 


number of starting element in a contour plot. 


E2 


number of finishing element in a contour plot. 


PRESS 


Same file as in MAPI. 


Notes on CONTI 





1. The X (axial) and Y (circumferential) dimensions of P must be entered 
before program is run. P[8, 128] is nearly the maximum array size that 
can be stored in the HP9830A. 

2. The axes are drawn so that X (axial) runs from K6 to K13 and Y (circum- 
ferential) runs from the start to the end of a blade pair. Contours 
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Figure 14. Flow diagram for CONTI 
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Table 8. Listing of CONTI 



j 0 R £ Pi WALL PREP SURE fl R F PROGRAM “CONI' i ’ 



20 L 


MIT 0 










30 r - 


ILES 


“RESS 








40 P 


EM IN 


3 EEC HI MENS IONS FOR P 


»:; t , j :• 


58 DIM PC 


3 3 1 £8 I 


* f 1 [ :< 3 ] ■ B C • J h o .1 !» . 


Zi 3 ]n QL 3 


1 0 ••• 


[ 1 ] = 0 










7 0 >: 


C 2 ] = u 


c- 
>; J 








S 0 '••• 


t 3 : * 1 










9 9 F 


C 4 0=1 


. -| -j 

■ z* r 








1 00 


ME 5 3= 


1 . 75 








1 l 0 


mi: s ]= 


2.1 £ 








120 


ME 7 3 = 


2. 5 








100 


Kt 8 1--- 


3 








140 


MC 9 3= 


7* cr 








150 


D 1 SP 


"ENTER 




! ! H X I fi L .•* fi H ]J 


V D I MENS 


1 6 y 


INPUT 


1 0 < . J 0 








1 7 0 


SCALE 


y h X [ y 


3 


a y B 5 j, J y t y „ 5 




.!. 3 0 


MAM I S 


0 . 5 3 4 


1 


0:24 [ 8 ] 




L 99 


MAM I S 


J 0 + 0 . 


5 


■» ? 0 i X [ 8 ] 




2 0 0 


7 fix IS 


MC 1 3< 


i 


■J y ° y . 5 1 J y / y . 5 




£19 


YAM I S 


MC 3 ] i 


i 


y i y n 5 ? 3 y y s !j 




220 


FOR I 


" 2 T 0 


•— i 
»' 






v>;m 


i fib I O 


MC I j 3 


i 


30 * 0 . 5 * ..J0 f y . 5 




240 


PEN 










250 


NEXT 


[ 








2 6 0 


READ 


# 1 3 1 








270 


M ri T k 


L H D tt 


i 


J P 




23Q 


LI 1 S r' 


ENTER 




CONTOUR P=" 5 




2 9 0 


I MF'IJT 


IJ 








3 0 0 


D I SP 


"ENTER 




FIRST & LfiST E 


LEi’lEN T 4 " 


3 1 U 


I NF'UT 


t. 1 3 E. £ 








320 


FOR E 


-£ 1 TO 




2 +* ' I 0 — I ) V J 0 - 1 


> 


2! 30 


i P E / 


E2 THE 


n 


1250 




343 


1 F E . • 


I NT i E/ 




y42 THEN 430 




350 


AC i i 1 


]= 1 NT ■ 


( 


E/2-0 - 0000 1 > 


J0-1 . ri 


3 b 0 


rlC 1 J £ 


J = E £~ 


C 


fit t « n-i >*< j0- 


i +■ 1 


373 


A C £ i 1 


J -AC i 3 


L 


.1 +■ i 




3 3 M 


AC £ j 2 


j = AC 1 3 


V* 


.3-1 




390 


AC 3 > 1 


I- AC 1 3 


1 


] +■ 1 




4 0 3 


AC 3 . £ 


]=AC 1 i 


£ 


.j 




410 


P ~ O 










423 


GOTO 


500 








430 


AC 1 ■■ 1 


> I N T •: 


v 


' E 4 1 ) •••*' 2 - 0 0 0 0 0 


L :* : jt; • ; 


4 4 3 


AC 1 j £ 


3=<.E-+ i 




2- fiC 1. fi 1 1 ' 


< jy - 1 > 


4 5 0 


AC £• 1 


J--AC 1 3 


l 


1 




4 6 0 


AC 2j 2 


3=flC t ’ 


£ 


3 + 1 




478 


AC 3 < 1 


1-HL 1 3 


1 


1 i- l 




430 


A 1 J 3 £ 


1 - A L 1 3 




•1 

.1 




4 9 3 


P 1 










530 


AC l - 3 


1-1 








5 1 3 


AC £ ’ 3 


.1= i 








520 


AC 3 3 3 


1=1 








5 30 


.:i i j= 


PC AC 1 - 


!. 


1 1 , : i .! 




5 4 y 


zi: £ i- 


PC AC £ ■ 


! 

. 1 , 


J* ML" Z £ ;i ] 
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•::c 3 i = 


r i. Ill 3 


1 


! 3 Hi . £ J 1 




5 b 0 


REM C 


t-lECi. i 


\ 


" 1 . :: ' I N i:: 




573 


if zi.: 


.1 j • 


f 


ri N O i . ■ i 


I zi!:: i ' 


5G3 


if ::c 


j: ] •■=• 


! ' 


and 3 ■ Z •. J 


T-i .i’ 


5 ii 


IF zc 


1 i - 




AND T -■ ZC '] 


1 ■! r ! • 


4 T i M 


j r , j 






a mi.) i; ■ i 


Ti 



I'iH : 



:■ Ul. 



hl 9 
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Table 8, 



Cont 



61.0 GOTO 1230 
620 MhT B= I HV (' H '■ 

630 MAT Q=B*2 

6 4 ti IF F’ v. 0i .b f HEN '? 4 ti 

650 IF QC 1 ] = 0 THEN 72 O 

660 X=< - QC 2 3*fll 1 . 2 3-QL 3 3 + C > /QC i 3 

670 Y=HC 1 j 2 3 

680 IF fl[ tstl <-- X HMD HC 3* 1 3 >* 4 THEN 
600 GOTO 720 
70O W= I MT ( X> 

7 10 PL 0 T X C W 3 + < X CW+1 3- X I M 3 > * ( X - W ) , V 
720 IF QC 2 3-0 THEM 700 
730 Y=< ~QC 1 3 + fiC 1 » 1 3 -QC 3 3-i-C>/Q[ 2 3 
740 X=fiC 1 » 1 ] 

750 IF fit 1 » 2 3 <‘= Y fiH Ii fit 2? 2 3 >= V THEM 
760 GOTO 700 
770 H= I NT 0-0 

780 PLOT XC M 1 H c. XC H+ i J-XC W .1 1 * r X ~N > , V 

? 0 U D C 1 !' 1 ] “ 1 

800 DC 1 ' 2 3= 1 

810 DC 2 s 1 j-QC 1 3 

820 DC 2’ 2 ] = QC 2 3 

830 FC 1 3 = fiC 3 j 1 1+fiC 3 1 2 3 

840 FC 2 J = C-QC 3 3 

850 M = DET C D 

860 IF M = O THEM 1230 

870 MhT G= I MV < D :> 

880 MHT H--G+-F 

800 IF HC 1 3 <= fit 3, l 1 HMD HC 1 3 = HI 2. ? 1 



000 


GOTO i 


230 








010 


M- 1 Ml ■ 


HC 1 3 :< 








020 


PLOT X 


C M 3 + ■ ; C 


i j - m i 


M 3 ;• -( , HC 1 3 1 1 


, 


030 


GOTO 1 


230 








040 


IF QC 1 3=0 I HE M 


1 0 i 0 






050 


X = ( -Ql 


2 J*fiC l > 2 3 


-QL 3 :i 


-C > Ql. 1 1 




0 6 0 


i -fit 1 ? 


2 3 








0 ? 0 


IF fill 


i 1 . ] •: X 


H * D H 


r 3,1 1 > - v '!' 


H h i' i 


080 


GOTO 1 


0 l 0 








0 0 0 


H= I MT '■ 


x .• 








1 000 


PLO C 


XL N 3 + - :■! M 


i i j-x 


L M me X-H ■ i' 




1 0 1 0 


IF QC 


2 3=0 THEM 


1030 






1 020 


Y = ‘" -Q 


I 1 1 -Hil. 2 ’ 1 


J 3 


J h: • Q[ 2 3 




1 030 


X--RC 2 


j 1 3 








1040 


IF RC 


2 1 2 J = V 


HMD 


HC3-2J >-- , 


The: 


1 050 


GOTO 


1 030 








1 060 


W= I II T 


'• r i 1 








:!. 970 


PL.O T 


1 t m i • ■ . ;r ' i 


t 1 J-ni 


!. M 3 ' * . X--H ' • 




1 080 


PC l > l 


]•-- 1 








1 000 


DC l , 2 


]=• i 








l 1 00 


DC 2 ’ 1 


‘J - Q C i 1 








1 1 1 0 


DC 2 1 2 


3 - Ql 2 J 









.1 i 20 FT I T' Hi 2 - , 3 -h 2> 2 
1 130 FC 2 -QC 'j ! 

1140 r"- DE • D 

1150 IF M-0 THEM J2 H 

1160 MHT [. • 

1 7 H MH T H--GF 
1 1 88 IF HC 1 I .'it. i ■ 1 

Li 08 i ,n ■ 18 y 1 
: .011 U III' Hi ! i 
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Table 8. Cont. 



1 2 1 6 


PLOT 


NE M ] + ( XC W+ 1 J--XH W 1 ) * HE 1 ]- U > i HE 


1220 


GOTO 


1230 


1230 


PEN 




1.240 


NEXT 


E 


125 0 


DISP 


" ANOTHER CONTOUR? l-YES* 0-NO" 


1260 


INPUT 


cs 


1 2 r 0 


IF 03 


>0 . 5 THEN £30 


1230 


END 


- 
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not quite go to these extremities because the pacer system starts 
sampling at 1/256 th of a blade pair and finishes sampling at 255/256th 
of a blade pair. Recall that the start of a blade pair need not coin- 
cide with a blade because of phase lags in the pacer system and the 
location of the timing disk relative to the rotor. Grid lines parallel 
to the Y-axis are the lines scanned by the Kulite transducers. 

3. If other transducer locations are used their axial distance downstream 
of K6 must be entered before the program is run. 

4. The variable plotted in the Y-direction is the column number in P. The 
variable plotted in the X-direction is the axial location of a transducer 
(downstream of K6) and is derived from the corresponding row number in P. 
When setting up the X-Y plotter it is advisable to make both the X and Y 
scales equal to twice full scale. Note that the circumferential distance 
along the wall across two blade passages is 3.847 inches. 

5. Entering of the first and last element numbers enables faster plotting 
of a contour which is known in advance to cover only a limited part of 
the field. 

6 . Triangular element representation of the surface defining the pressure 
coefficient distribution, Cp (X,Y) 

The surface (X,Y) is approximated by triangular elements as 

illustrated in figure 15. Element numbers E are as shown. This process 

represents linear interpolation between measured pressure coefficients. 

Contours are obtained from the intersection of planes C = constant with 

P 

this approximation to the pressure distribution Thus the contours are 
composed of straight line segments. The nodes correspond in the Y-direction 
(circumferential) to points at which pacer data is available. The node 
numbering, in the local sense, for typical odd and even numbered elements 
is shown in figure 15. 
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Figure L5 Numbering of triangular elements in CONT and CONTI. Also 

numbering of nodes (vertices of triangles) for odd and even 
numbered elements. 
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If Z = aX + bY + c is the plane containing the three nodes of a 



triangle then 



z(l)" 




n yi l 




r *+ 

a 


Z(2) 


= 


X2 Y2 1 




b 


_Z(31 




X3 Y3 X. 




.C- 



where Z(I) is the pressure coefficient at local node number I (i.e. at 
XI, YI). 

This matrix equation is written as 

Z = Aa 

and hence 

Q = A = BZ say. 

If E is odd , the node coordinates are related to the element number 
as follows: 

XI = A[l, 1] = INT ( ( (E + l)/2 - .00001) / (JO - 1) ) + 1 

Yl = (E + l)/2 - (A[l, 1] -1) * (JO -1) 

X2 = XI , Y2 = Yl + 1 

X3 = XI + 1 , Y3 = Yl 

Note that the number .00001 is included to avoid problems associated with 
round-off error. 

Z = aX + bY + c is plane containing triangle. For Z = c', line of 
intersection with triangle (i.e. contour) is given by aX + bY + (c -c') = 0. 
On the side X = I = A[l, 1], (See Figure 15) 

Y = (- aA[l, 1] - (c - c'))/b 

= (- Q [1] * A[l, 1] - Q(3) + c')/Q(2) 
provided Q (2) # 0 
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If Q ( 2) = 0 , the contour is parallel to the side in question. On the 

side Y = J = A[l, 2] 

X = (-Q[2] * A[l, 2] - Q ( 3) + C f )/Q(l) 
provided Q (1) £ 0 

On the hypotenuse 

X + Y = A[3, 1] + A[3, 2] 
aX + bY = c 1 - c 
Hence 



1 1 


f X l = |a[3, 1] + A[3, 2] 


a b. 


i — 

0 

1 

o 

1 



This matrix equation is written as 
DH = F 

Hence H = GF where G = D ^ 

l-Srf <*'*<' 

If E is even , the node coordinates are related to the element number 
as follows: 

XI = A[l, 1] = INT ((E/2 - .00001)/(J0 - 1) ) + 1 

Y1 = A[l, 2] = E/2 - (A[l, 1] - 1) * (JO - 1) + 1 

X2 = XI + 1 , Y2 « Y1 - 1 

X3 = XI + 1 , Y3 = Y1 

Z = aX + bY + c is plane containing the triangle. For Z = c 1 , line of 

intersection (i.e. contour) is given by aX 4- bY + (c - c 1 ) = 0 . 

On the side X = I = A[2, 1] , 

Y = (- Q(l) * A[2, 1] - Q(3) + c’)/Q(2) 
provided Q ( 2) # 0 

On the side Y = J = A[l, 2] 

X = (- Q(2) * A [1, 2] - Q(3) + c’)/Q(l) 

provided Q(l) j 4 0 . 
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On the hypotenuse 

X + Y = A[2, 1] + A[2, 2] 
aX + bY = c' - c 



Hence 




A[2 , 1] + A [2 , 2] 



a 




- c 



This matrix equation is written as 
DH = F 

Hence H = GF where G = D ^ 

A/ ^ 



The intersections of the contour with the sides of the triangular 
elements are computed in the manner outlined above. Tests are performed 
to identify those points of intersection which lie on side of the triangle 
as opposed to intersections which lie on extrapolations of the sides. 
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APPENDIX H 



DETAILS OF CONT 

CONT is used to plot contours of constant casing wall pressure from 
an array of pressure coefficients in which both the axial and the circum- 
ferential spacings are uniform. It can be used to plot contours for an 
array output by MAP2. The program differs from CONTI only in its plotting 
of x-coordinates of contours. A listing is in Table 9. 

Notes on CONT 

1. If MAP2 is used for interpolation, the file PRESS must be replaced by 
file INTER and P must have dimensions P[29, 32] . Contours must be 
plotted by processing the array of interpolated pressure coefficients 
in four sections. This requires shifting of origin on the x-y plotter 
for each quarter. 

2. The variables plotted in the x and y directions are respectively 

the row number and column number of the sub-array of pressure coefficients. 
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Table 9. Listing of CONT 



1 

». r -' 


REN WALL PRES 


SURE NAP PROGRAM "CONT" 


15 


UNIT 0 




IS 


FILES INTER 




20 


PEN #30 DIM P 


Ifh 10 ) 


30 


DIM PC 25 > 32 3 > 


H C 5 ?» 3 .1 ■» B [. 3 ? 3 j * Z [ 3 ] ? Q C 3 1 -» r - 


4Q 


CUSP "ENTER X 


AND V DIMENSIONS" 5 


50 


INPUT 10- J0 




60 


3 C RLE 1 ' I U ' U ? 


J0 


70 


x rx is 0 m: [ 0-1 


• 1 0 > 1 J I 0 


8 0 


XflXIS J0J (10- 


1 > 10- 1 - I 0 


90 


T FIX I S' 1- 100-0 


4 J0 


1 00 


YfiXIS 10- 100 


;i 0 ;i J 0 


110 


READ #1 > 1 




120 


MAT READ # 1 


! P 


170 


Ii ISP "ENTER 


1 ONTOiJF. P~" ; 


ISO 


INPUT C 




1 35 


DISP "ENTER 


FIRST ft: LAST ELEMENT#"" 


1 O O 


INPUT E 1 - E2 




1 90 


FOP E-El TO 


2 * ( 1 0 — i » •+• ( ... 1 0 •• 1 j 


1 9 5 


IF E >E2 THEN 


i 0 U € 


2 8 0 


IF E> INK E/2 


..•*2 THEN 290 


2 1 0 


AC 1 f 1 KINK ■ 


E 2 -- 0 „ 0000 1 .< ( J 0 l r [ 


220 


AC t - 2 ]=E 2-'. 


AC [ 1 i J-i '<*•; JO - 1 :• .1 


2 30 


AC 2* 1 ] = AC 1 - 1 


3 + 1 


240 


H [ 2 9 c ] : = H !. 1 ? 2 


1-1 


250 


AC 3 « 1 ] = flC 1 j 1 


.1 1 i 


2 6 0 


AC 3 1 2 ]=rlt 1 ■> 2 


j 


270 


R - 0 




230 


GOTO 360 




2 9 0 


AC i - I KINK < 


E r ! .' 2 - U i t-0000 i 1 . 0 - ' • '■ 


3 0 0 


AC 1 - 2 ]■=( .EH : 


2 -'-.AC 1 j 1 3-1 ’ M..JU- 1 ■ 


3 1 0 


A L d ? 1 .1 - A C 1 ? 1 


j 


320 


AC 2 - 2 KAC 1 - 2 


3-i 1 


3 3 0 


H C 3 ' i K A L 1 ’ 1 


3+1 


3 4 0 


H C 3 - 2 K H C l ’ 2 


1 


350 


R= 1 




3 6 0 


ACK 3 K l 




• " “7 i T-i 
i r.» 


AC 2 - 3 ]- 1 




3 8 0 


AC 3- 3 1-1 . 




h* y y 


ZC 1 3 -PC AC 1 - : 


' + r i < 2 j j 


J 1 |7i 
*-r x U 


2C 2 3= PC AC 2 - 1 


] « HI. 2 * 2 J J 


4 2 0 


2C 3 3 - PC AC 3 - 1 


m h[ 3, 2 3 1 


421 


PEN C. NEC 1 ' IF 


C IS IN E 


422 


I F Zf 1 3 •• -- c 


AND 0 - Zf ... ] T-i,- ri 4 - 


4 2 3 


IF zr 2 J - 1 . 


AND C /= 2T i. j 1 HEN 


4 24 


IF zr 5 j ... 


AND !" •= EC I- .! TAR F : . ' 


425 


:. F Zl. 3 ■ .. 


-i'i H L •- EC L 3 " HEN ! ’ 


h *' !*s 


30 ‘0 ! 0O A 




..j. 


NAT A IN V • - ■ 




4 > y 


MAC OB*Z 




4 4 0 


• .• - R i . r ~i 1 1 j : 


1 • I !■' 


4 4 5 


-F Oi !. 3---M i'i 


l:: 1 i »! 


4 5 0 


:k- -in 2 3 - ri ; 1 


. ■ :■ -! 3. i ~ 1 «. • i. ■ J 


4 6 0 


"-hC i - 2 3 




4 70 


• f R r j . i i 


IN j AC in i ■ 


•h- !"■ 0 


• TO '>10 




•7 : .n' i 


i- • t r v • i 





"j - I. ! ; M 1 i 






Table 9. Cont. 



505 Y-' -QC 1 3*AC i ? 1 3-QC 3 3+C)/QC 2 ] 

510 K=--RC 1-13 

520 IF AC 1*23 <= Y AMD At 2 .- 2 3 >= Y THEN 

530 GOTO 57Q 

540 PLOT X < Y" 1 

570 DC 1 1 1 ]=1 

580 DC 1 < 2 3=1 

590 DC 2 ' 1 J=QC 1 3 

680 DC 2? 2 3=QC 2 3 

610 FC 1 3=AC3i 1 3 + fit 3*23 
620 FC 2 3=C-QC 3 3 

625 M=DET ( D > 

626 IF M=0 THEN 1088 
638 MAT G= I HV C D > 

6 4 8 MAT H=ij*F 

645 IF HC 1 3 < = fit 3-13 AND HC 1 3 >= AC2) 1 3 

646 GOTO 1080 

638 PLOT HC 1 3 j HC 2 3- 1 

660 GOTO 1800 

678 IF QC 1 3=8 THEN 758 

6 9 O X = < -QC2 3*fiC 1)23 - Q C 3 3 + C > / Q C 1 3 

700 Y=A[ 1 ) 2 3 

718 IF HC 1)13 < = X AND AC 3)1 3 >= X THEN 

720 GOTO 750 

? 3 0 P L 0 T X j Y - 1 

758 IF QC 2 3 = 0 THEN 838 

768 Y= C -QC 1 3* AC 2 j 1 3-QC 3 3+C > /QC 2 3 

770 X = AC 2)13 

788 IF fit 2) 2 3 Y AND fiC 3)23 >= ■' THEN 
798 GOTO 830 
880 PLOT M'Y”! 

838 DC 1) 13 = 1 

848 DC 1)2 3=1 

850 DC 2) 1 3 = QC 1 3 

3 6 8 DC 2 1 2 3 = Q C 2 3 

870 FC 1 3=fiC 2) 1 3 + AC 2)2 3 

836 FC 2 3 = C-QC 3 3 

885 fl = D E T 1 D > 

836 IF M=8 THEN 1088 
890 MAT G = I NV C D > 

908 MAT H=G*F ■ 

985 IF HC l 3 >= fiC 1)13 AND HC 1 3 •• = HC2? 1 3 

966 GOTO 1008 

9.10 PLOT HC 1 3 1 HC 2 3- I 

920 GOTO 1000 

1 0 O O P E N 

.1005 NEXT E 

1 0 0 6 D I S F' " A N 0 7 H L P C 0 N T 0 i j k ? 1 • • Y E ■ 0 -• ■ 0 ’ 
1007 INPUT C 3 " 

1 0 8 8 I F C 8 > 0 .. 5 T H E N i 7 0 
1810 END 
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APPENDIX I 



DETAILS OF PLOTS A AND PLOTS B 

PLOTS A and PLOTS B are used to plot pacer raw data and free-run 
raw data respectively against circumferential distance • The programs are 
listed in Tables 10 and 11. Note that PLOTS A, with I = 1 in line 110, 
plots the output of the designated Kulite transducer. If I = 4 the raw 
one per blade signal on A/D converter channel number 0 is plotted. 
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Table 10. Listing of PLOTSA 



10 


REM — 





* •+• $ 


-Kp 


L 


| j “j“ Q |Zj -A 


20 


UNIT 


0 










30 


FILES 


REC 


v # ? 


DA 


Y 


AVI 


40 


DIM D 


SC 5 j 


1 28 


] 






50 


DISF 


“ENTER 


r L 




#" ; 


b 0 


I NF'IJT 


N 










70 


RE FlD 


#2 j N 










80 


MAT R 


EAD 


4 2 


; d 






80 


SCALE 


P < 


M0 ' 


- 1 




£T -j cr 


1 00 


nisp 


"DR 


AN 


h-: 


E 


s? ^ ’ i =VES 


1 10 


INF'IJ 


T NO 










.1 20 


IF N 


0 = 0 


THE 


N 


I 


50 


1 30 


XflX I 


S @ J 


5 O i 


0 .i 


3 


00 


140 


VAX I 


s 0 , 


0. 3 


i "■ 


i 


.5*1.5 


1 j u 


1 = 1 












1 60 


FOR . 


J= l 


TO 


12 


3 




170 


PLOT 


2* J 


■■ DC 


I j 


J 


] 


180 


NEXT 


...I 










190 


PEN 












200 


END 
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Table 11. Listing of PLOTSB 



LT 


Erl " 


PLOTSB " -$ 


If-! MOMS- - 


1 @ 


UHI T 


0 




'cl 0 


FILES 


RECY# , HR 


TRY! 


3 0 


DIM D 


SC 16’ 102 ] 




40 


D I SP 


"ENTER RE 


L li ' ? 


50 


INPUT 


N 




h 0 


REfiD 


#2 1 N 






MAT R 


EAD # 25 D 




30. 


SCALE 


0 1 1 ?00 , - 


1 . 5 1 1 . 5 


31 


n i sp 


" DRAW AXE 


bY 1=YE 


I' ul 


INPUT 


NO 




33 


IF M0 


=0 THEN 1 


10 


?0 


NAN I S 


0 1 00 ’ 0 ’ 


1700 


100 


VAN I 


S 0 ’0.5?— 


1 - 5 * 1 . 5 


1 1 0 


FOR 


1 = 1 r o 1 6 




120 


r U R 


J=1 TO 10 


1 


1 3 0 


PLOT 


r T — 1 -f [ u 


1 i j. m t , 


140 


NEXT 


J 




15U 


NEXT 


i 




155 


PEN 






1 69 


END 
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APPENDIX J 



DETAILS OF TITIPK 

TITIPK is used to superimpose blade tip profiles on wall pressure 
contours. The program is loaded onto the programmable keys of the HP9830A. 
The program is listed in Table 12. The blade tip profile is tabulated by 
PAIGE (1976), figure 2 and table 1). The axial and circumferential units 
are inches and the program is compatible with a wall pressure map which has 
dimensions of 3 inches axially by 3.847 inches circumferentially. The 
"lower left" and "upper right" on the plotter should be set to the corre- 
sponding points on the wall pressure map. 

The key programs should be "continued" after CONTI has been run so 
that the Kulite data is available in main memory. When <CONT> <f^> 
issued, the location (YO) of the blade leading edge from the lower boundary 
of the wall maps is calculated from the data of transducer K9 and appears in 
the display. <CONT> <f^> is then issued and YO is requested as an input. 

The blade profiles are then drawn and the key program ends. The contour 
plotting with CONTI can be continued by issuing <CONT> 1250. 
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Table 12. Listing of T1T1PK 



i. .1 Ki. .11 



• 1 i : 


H! ■ 


• i" l: ■/ 


pr u 


Of! Ml I 


PL 


;.j II 


!. --0 










‘10 


i./hii; ■ ! 


- : i u 


• J , • . ; 4 






50 D 


RT'h .. ' 


T i !;;;! A :* 


t 6 7 j 


\ F.i B '0 i i ! :..i s ! 


j 


63 F: 


END i'i:i 


- RO ,i 


! 1 1 i !• ! .. .! 


• ! ! 0 




70 U 


7 3P ‘ i ! 


•sir R 


’r'i, " 






«*:i i 


hPu r V 


r i 








V 


0-Y9 










j. 00 


I) EG 










1 1 0 


R 1 = »: l n 


- F’ 


1 * F 


0 F' 1 




J <20 


ru=fi r i-i 


* li 1 


3 i-i 7: 1 


-Hi 42 




1 00 


ji i -a i •- 


G 1 - 90 






1 4 0 


V 3 0 










1 L j0 


V.3-R i 










1 60 


;•••: 4 =xs+ 


•: c 0 •• 


2HM - 


rC OSG 1 




170 


4-Y J + 


1 C0- 


2-J R i ) 


4SIHG1 




100 


! 1 5 « ';■■■, j 


. R 0 - 


R 1 ’■ * C 


OSD 1 




1 ,; ^0 


V5~Y 3 i 


■ PO - 


R 1 •) if S 


1 NT 1 




2 0 fi 


REM - - - 


- - F F : 


0 F ILL 


T I P ( NB 


■i - 


: j 0 


1 3 ■“* 6 










220 


r 3-90 - 


jji 








2 30 


T 4 = ■: T 3 


T 2 > 


£ kl 






240 


FOP T 1 


r 

• t c.» 


TO 12 


•3TEP T 


4 


2 


3 • p 


1 4 q r 


NT , 






2 60 


V - Y 3 R 


1. *C0 


3T 1 






.t 7 0 


IF Y + 'i 


0 • 3 „ 


34 7 r 


HEN 310 




2 ft hi 


L F i * 


0 0 


THFN 


3 1 0 




20U 


^ LOT 


t ; ; 0 * 


Y + Y0 






3 00 


! jQT ij ; 


2 0 








3 1 0 


PEN 










3 2 0 


NEXT 7 


1 








3 3 u 


PEN- 


- ••• 3 ij 


6 T lull 


! -IDE- 6 


c > 


340 


T2=-U:! 










350 


r 3-2 !" i 


i - D 1 








3 60 


T4=i Tj 


- r 2 V 


• 200 






3 7 0 


f'OR T !. 


•- r 2 


7 0 T7 


■ STEF 14 


380 


X-V5 -ptj t- f:ij 


ST 1 






39G 


• - Y5 + R 


0-3 I NT l 






400 


lF V + V 


0 3 3 


34 7 1 


HEN 44ti 




4 j. 0 


[F ,'i-V 


0 ' 0 


THEN 


440 




420 


plot 


i- ' ! r j ;i 


V +■ V 0 






430 


!.iU f 0 4 


5 0 








440 


FEN 










450 


NEXT ; 


i 








460 


REM— - 


— P F 


OF ; Li 


T . E . ! L 


Ti •: 

4.' • 


4 7 0 


i i J ~ ~ l.j v 


+ R I 








480 


T 3 = I 3 6 


i-Gi 








4 9 0 


T4 = T 3 


i - T 2 


.• 20 






5 S 0 


FOR n 


- T 2 


TO r; 


■ 3 T E P T 


’ -* 

■t 


5 i 0 


:-=:=x 4 


: 1 * 3 1 


N r , 






5 2 0 


Y = 'Y 4 t- F 


' l C C 


i:3Tl 






530 


i F Y + 


' 0 3 . 


347' '! 


"HEN 570 


1 


540 


j. i i r i 


' 0 < 0 


THEN 


570 




550 


p * L r 


; + ‘ ; ij 4 


V 4 V 0 






5 b 4 


GO TO 2 


1 8 k;i 








570 


PEN 










580 


NEXT i 


1 








59 U 


21 1 — X 











6 0 0 ' ' 1 = V 




1 ■’ 



4 
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Table 12. Cont 



t i L 


; 1 1 - 


pi. 


eg ..if- j. u r . d 


r , ’ : 3 


I j; =- 


‘ E L •• 


; iiv. < 


P 3o 


Y2~Y7 


-R 1 1 ' 


T S G : 


r ■4i-;1 


FOR I 


1 ! P 


i V.' ’ 


b 3 0 


N~N 1 


1* 1 

i i. •' 


r •' , ' 1 j ... ROO 


b b 0 


i j ^ ~r 


‘ N 32 


'■ +• 'T 1 - v e v :•••: i 


670 


1 F t‘ "f 


Y0 -3, 


04 , THEM 71 a 


6 o y 


IF Y 4 


Y0< 0 


THEN 710 


6 9 0 


PLOT 


■ •, , 0 h 


■; + 'V'0 


700 


GOTO 


7 2 0 




7 1 0 


PEN 






7’ 0 0 


NEXT 


!' 




730 


PEN 






7 40 


IF N , 


:• 1 . 5 


THEM 706 


7 3 0 


N 1 “H 1 + t 




T 6 !j 


\ 0 - Y 9 


1 , >2 


0 5 S-2*N1 > 


770 


laO TO 


200 




780 


END 







10 


-TH-- 


- 


-POUT I ME 




0 - EG 8 4 


7 




43 


FOP 1 = 


i. 


TO i 28 


ri 0 


IF PC 4 




I K 8 THE 


h", 0 


V 1 -PC 4 


* 


I J 


7 0 


GOTO 1 


4 ij 


£ 0 


IF 1=1 




THEM 1 -1-8 


90 


I F V 1 


0 


THEM .14 


1 3 0 


V 2= PC 


■) 


’ i 1 


1 1 0 


i 0~ H ■ l l 


( 


1-1= 5tV1 


1 23 


r« i sp 


■ ' 


T8-- "T8 


1 30 


STOP 






40 


NEXT 


i 




1 5 0 


END 







EDGE FOR T I T IRK 



-on y 
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APPENDIX K 



OPERATING PROCEDURE FOR DATA ACQUISITION SYSTEM 

1. Load Real Time Executive Basic into HP21MX computer. 

2. Load KULITE into HP21MX. Tune pacer. 

3. Put disk labeled "Transonic Compressor - Paige" into HP9867B mass memory. 

4. At start of a run, get RESET1 from unit 0 and run RESET1 to initialize 
the number in file RECY// to unity. 

5. Scratch RESET1. 

6. Get TRAN4 from unit 0. 

7. Run TRAN4. The HP9830A display will remain blank while the HP9830A 
awaits data from the HP21MX. 

8. Run KULITE, noting that the two mode switches on the Pacer must be set 
according to the A/D converter mode to be used (i.e. 0 or 4). In cali- 
bration, four scans of the twelve channels must be made. The first scan 
must be with the pressure at S2 on the reference side of the Kulite 
transducers. The second and third scans are made with other steady 
calibration pressures on the reference side. The fourth scans must be 
made but any signals can be used on the A/D converter provided that they 
do not cause overloads. In mode 0 the pacer must not be altered during 
an experiment. 

9. On completion of an experiment both the HP9830A and the HP21MX must be 
reset as instructed by their displays prior to performing another 
experiment. 

10. On completion of a run the data which is stored temporarily in DATAY1 
file must be transferred to permanent files CKRWra duplicated on both 
unit 0 and unit 1. First open CRKWm where m is the run number. The 
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length of CKRWm must be set at k records where 
k = record number for last experiment + g 
and g = 1 if last experiment is a calibration 

5 if last experiment is a pacer experiment 

13 if last experiment is a free-run experiment 

11. To abort the HP21MX program, enter AB. 
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